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A growing or shrinking disc will adopt a conical shape, its intrinsic geometry characterized by a 
surplus angle ipe at the apex. If growth is slow, the cone will find its equilibrium. Whereas this is 
trivial if ipe < 0, the disc can fold into one of a discrete infinite number of states if ipe > 0. We 
construct these states in the regime where bending dominates, determine their energies and how 
stress is distributed in them. For each state a critical value of ipe is identified beyond which the cone 
touches itself. Before this occurs, all states are stable; the ground state has two-fold symmetry. 

PACS numbers: 68.55.-a,46.32.+x,02.40.Hw 
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Soft matter systems may display enormous complexity at 
the microscopic level [l|, |2|. However, on mesoscopic or 
larger scales of physical interest, the relevant degrees of 
freedom very often turn out to be purely geometrical. If 
one dimension is much smaller than the two others, an 
effective description in terms of a two-dimensional surface 
becomes appropriate. This is just as true for biological 
membranes [3, 4^ as it is for inanimate matter [5] . 

The equilibrium shape of the surface is often a minimum 
of bending energy. Typically however, one must take into 
account external forces, or constraints on the geometry: 
these may be global as in the fluid membranes occurring 
in cells in which the area or the enclosed volume is fixed; 
they may also be local as in plant tissues which are de- 
scribed, to a good approximation, as an unstretchable 
surface with a fixed metric [6|, [D] . 

In general, bending is not possible without stretching. It 
turns out, however, that the most effective way to mini- 
mize the energy and satisfy the constraint is by confining 
the regions where stretching occurs to a series of sharp 
peaks and ridges |8','9']. On length scales much larger than 
the thickness of the sheet, these can then be treated as 
points and curves along which boundary conditions are 
set on the surface. The simplest geometry of this kind is 
the developable cone-the "point defect" of folding, with 
the bending energy localized near the apex. Such a ge- 
ometry is illustrated beautifully by a flat planar disc of 
paper depressed into a circular frame by applying a point 
force to its center [10|. 

Conical shapes also occur in living tissues. The unicel- 
lular algae Acetabularia acetabulum, for example ,_grows 
a conical cap in the course of its development [71, UjJ. 
Despite the superficial similarity, however, the point-like 
singularity exhibited in such cones is very different from 
the one which appears in the developable cone. In the 
latter, the singularity at the apex is extrinsic; the met- 
ric itself remains the same as that of the original disc. 



This is captured by the fact that the surplus angle at the 
apex vanishes: the length of the closed curve at a unit 
distance from the apex is equal to 27r. If however there 
is a surplus or a deficit, there will be a non-trivial folded 
state even when external forces do not act. Whereas this 
state is an unremarkable circular cone in the case of a 
deficit, when the deficit is turned to surplus, the folded 
shape-an excess-cone (e-cone for short)-exhibits a sur- 
prisingly subtle behavior. In this letter, we will describe 
the equilibrium states associated with these "point de- 
fects" in the full nonlinear theory. Remarkably, there 
exists a discrete infinity of states (they are "quantized"); 
they can therefore be completely classified in terms of 
the surplus angle ipe and a quantum number n. For each 
n we show that there is a critical value of (fe, increas- 
ing monotonically with n, beyond which the cone makes 
contact with itself; if cpe > 35.23 self-contact cannot be 
avoided. The resulting shapes are not unlike the collars- 
called ruffs-one associates with portraits by Rembrandt 
or Hals. 

To understand the physics of e-cones we identify the 
stresses which underpin their geometry. We also address 
the question of stability. What is the ground state? Are 
all equilibria local minima of the bending energy? To an- 
swer these questions, one must take care not to lose sight 
of the local constraint of isometry. It is reassuring that 
one can easily built paper models and put the results of 
this analysis to the test. 

Geometry and energy. The e-cone can be parametrized 
in terms of a closed curve F : s ^ u{s) on the unit sphere 
where the arc- length s runs from to Sg = 27r + cpe along 
this curve (see Fig. [T]). If r denotes the radial distance 
from the origin, the surface is described by the vector 
function X(r, s) = ru{s). Its direction in Euclidean 
space R^ will be given by the polar and azimuthal an- 
gles on the sphere, (p and i^ respectively. The tangent 
vectors to the e-cone are u and t = u' where the prime 




FIG. 1: Geometry of the e-cone with ip^ — '^ 



denotes a derivative with respect to s. Together with the 
normal n = u xt these vectors form a right-handed sur- 
face basis. In the direction of u the surface is flat. The 
curvature along s is given hy tz = —n - t\ It is easy to 
show that hz is also the geodesic curvature of F on the 
unit sphere. It should, however, not be mistaken for the 
Frenet curvature. 

Our first task is to identify the configurations that mini- 
mize the bending energy of an unstretchable cone of ra- 
dius R with a surplus angle (pe > 0. If we introduce the 
cutoff ro at the apex and integrate over the radial direc- 
tion, the bending energy is given hy B = (a/2) ^p ds k?] 
the dimensional dependence of B is absorbed into the pa- 
rameter a = In {R/vq). While it will set the magnitude of 
the stresses in the cone, it does not play any direct role in 
determining the shape of the e-cone. The constraint of 
unstretchability is implemented by adding a term to the 
energy functional which fixes the metric via a set of local 
Lagrange multipliers T^^ [12|. These can be identified 
with a conserved tangential stress. 

The shape equation and its solution. There is a remark- 
ably simple way to determine the shape of the surface: 
first recall that to every continuous symmetry of a system 
a conserved Noether current exists. As the apex of the 
e-cone is fixed, translational invariance is broken. The 
bending energy is, however, rotationally invariant. The 
corresponding conserved vector J, related to the torque 
about the apex, is given by [l2| 



J /a = [k^I2 - C\\ 



Hi't- 



(1) 



We will suppose that J can be aligned with the z axis, 
J = Jz. Its square directly yields the first integral of the 
shape equation of the surface 



J' - c, 



+ KV4 + (i-q|)K2 



(2) 



where we define J := J /a. The constant C\\ is asso- 
ciated with the fixed arc-length; it will determine the 
stress established in the surface. Eqn. ([2]) is identical 
to the equation describing the behavior of planar Euler 



elastica with (scaled) tension a := C\\—l and n in place 
of the Frenet curvature. It is completely integrable in 
terms of elliptic functions [l3|, [ij] . The intrinsic closure 
condition ^.{sq) = nijS) = will provide a "quantization" 
of the solution: the e-cone has to be periodic in equi- 
librium with a period of 5e/n, where n is the number of 
folds. Solving Eqn. ([2|) for n we obtain 



k{s) 



-k [lC{k)/S] sii[2slC(k)/S,k\ 



(3) 



where S = Se/2n. The function sn (s, k) is the sine of 
the Jacobi amplitude am(5,/c) with parameter k. The 
symbol ]C(k) denotes the complete elliptic integral of the 
first kind [l5| . The parameter k is directly related to the 
stress C\\ and the torque J. 

To obtain the shape, consider the projections of J with 
respect to the local trihedron. Projecting onto u 



J{Z • U) = J COSl^ = K. 



(4) 



yields the polar angle i^ as a function of s since f<i{s) is 
known. Eqn. ^ places a strong constraint on the equi- 
librium shape. In particular, it implies that the only 
equilibrium shape consistent with a deficit angle is a cir- 
cular cone. Projecting J onto t reproduces the derivative 
of Eqn. (|4j). The remaining projection 



J{z -n) = Jlp' sin^ ^ = n^ /2 - C\\ 



(5) 



allows one to determine the azimuthal angle Lp via a sim- 
ple integration. 

The closure of the surface in Euclidean space R^ provides 
a second {extrinsic) closure condition, i. e., (^{sq) = 27r. 
This equation identifies the constant k implicitly as a 
function of the surplus angle ipQ and the quantum number 
n. Its numerical solution can be approximated remark- 
ably well by a polynomial of third order in (f^. For small 
(Pe, we find k = aiif^ + 0{ifl), where ai = -^(1 - 7^). 
Surface shapes. The e-cone has to have two or more 
folds. This is a consequence of the four- vertex theorem 
[l6|,ll7|. If the surplus angle is small, one finds a solution 
with no self-contact for all natural numbers n > 2. As 
an example, the first three n-folds for cpe = 27r are plot- 
ted in Fig.[2l^b)-(d). It is child's play to construct paper 
models; the 2-fold illustrated consists of two circular pa- 
per discs, each with a radial cut, glued together along 
the opposite sides of the cut. The model [see Fig. Efa)] 
closely resembles the calculated shape even though the 
surplus angle in question involves deformations of the 
flat geometry well outside the linear regime. 
As (fe is increased, the conical geometry becomes more 
crowded and, at some point, the mathematical surface 
will intersect itself. This happens first with the 2-fold. 
The physical surface, of course, will not self-intersect. 
Where different regions come into contact, they will ex- 
perience forces and they will deform accordingly [18]. To 
determine the critical surplus angle (p^^^^ above which 






FIG. 3: Stress C\\ as a function of (fe for n — 2 (solid line), 
3 (short-dashed), 4 (dash-dotted), 5 (long dashed), and 10 
(bold solid). The red curve shows C\\ for the touching 2- fold. 
Above C\\ = 1 the tension a along the curve V changes sign. 



FIG. 2: Paper model (a) and calculated surface shapes for 
Lpe = 27r with n = 2 (b), n = 3 (c), and n = A (d). 
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29.38 


34.92 


35.23 



TABLE I: Kissing points for different n-folds. 



this happens, consider the opening angle a = ^{S/2) at 
the turning point. It is given by a = arccos (J~^A>:max) 
with A^max = ^\f—kKL{k)/S. When a = the two- 
sides of the 2-fold touch along the z axis and Eqns. (|4|) 
and ([5]) simpHfy to J = /^max and J^ = 2C||. This 
implies that S'^ = 4(1 — k)JC{k)'^. Solving this equa- 
tion together with the extrinsic closure condition numer- 
ically yields k ^ —0.28 and (p^^^^ ^ 7.08. For e-cones 
with more than two folds, the situation is more compli- 
cated. Adjacent folds touch pair wise at some nonvan- 
ishing opening angle. The corresponding kissing condi- 
tions are now ^\sl.-^^J = and V^(s^iss) = (2j — l)7r/2n 
where j G {1, . . . , 2n}. In Tab. [J the values of (p^^^^ of 
various n-folds are given. Interestingly, cp^^^^ converges 
to (Pe^max ^ 35.23 from below if n is sent to infinity. 
This implies that one cannot find a stable surface with 
(fe > (/:'e,max which docs uot make contact with itself. 
The detailed analysis of these touching geometries lies 
beyond the scope of this letter. The 2-fold, however, is 
straightforward to study since it will touch itself only at 
two segments of the unit circle of arc-length (pt each. Us- 
ing the kissing conditions with n = 2^ (pt and C\\ can be 
determined simultaneously for any given (Pq. 
Stresses in the e-cone. Even though there are no exter- 
nal forces acting on the cone, stresses will be set up in 
the surface due to bending as well as the constraint on 
the metric. The stress tensor T^^ which fixes the latter is 
purely tangential. It is diagonal with respect to the basis 
(n, t) and constant along curves of constant radial dis- 
tance r. Its nonvanishing components along t and u are 



given by Ty = — Cy/r^ and, for a large disc, T± = — Ty 
[l2|; it is non-isotropic. 

In Fig.[3]C|| is plotted as a function ofcpe for different val- 
ues of n. For small values of cp^ the expansion k ^ aicpe 
can be used to write C\\ = (1 — n^) — (47r)~-^(3 — 7n^)(/:?e + 
0{(pl). In this regime the stress C\\ is negative; this cor- 
responds to a compressive stress along the tangential di- 
rection; an equal compensating tensile stress will act ra- 
dially. C|| is non- vanishing when (pe = 0^ representing 
the critical compression necessary to buckle the planar 
sheet into the corresponding mode. For a fixed surplus 
angle the absolute value of Cy increases with n. If (pe is 
increased, the curves for different n converge and appear 
to intersect in a single point (the 2-fold is exceptional 
making contact with itself before this point is reached). 
Investigated more carefully, however, a set of adjacent 
intersection points is found which converge to (Pq ^ 7.47 
for n ^ oo. Above this region each curve reaches a max- 
imum which diverges quadratically with n. If n > 5 and 
(pe sufficiently large, Cy is greater than 1. This implies 
that the tension a in the corresponding planar Euler elas- 
tica changes sign [of. text below Eqn. ([2])]. 
However, one must remember that the full stress in the 
e-cone includes a contribution due to bending. This 
becomes increasingly important as cp^ gets larger. The 
tangential projection /y of the full stress tensor can be 
written as /y = {J/r'^) sini^cp^ where cp is the basis vec- 
tor of (p. The transmitted force per length along F is a 
maximum at the equator and a minimum at the turn- 
ing points. It always points in the direction of cp; the 
tangential part of /y becomes tensile where (p' is nega- 
tive. These results can easily be verified by cutting the 
paper model along the flat direction and observing how 
the sheet reacts. 

Bending energy. We are also able to obtain an analytical 
expression for the bending energy B := B/{aSe/2) using 
the expression for k in terms of elliptic functions: B = 
64/C(A;) [S{k) - JC{k)] {n/s^f. We have normalized the 
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FIG. 4: Scaled bending energy B for n — 2 (solid line), 3 
(short-dashed), 4 (dash-dotted), 5 (long dashed), and 10 (bold 
solid). The red curve shows B for the touching 2- fold. 



bending energy by dividing by the area of the e-cone. 
The function f (/c) is the complete elliptic integral of the 
second kind [15]. 

For fixed (^^^ B scales essentially with v? since k depends 
only weakly on n. In Fig. 2] the energy is plotted as 
a function of Lp^ for the lowest n-folds. For small (/^e, 
all curves behave as 5 ?^ \^^ ~ l)v^e- The 2-fold is 
the ground state if (pe is below the critical value associ- 
ated with self-intersection. This observation is confirmed 
nicely by the paper model [see Fig. [21(a)]. For higher val- 
ues of (fe one needs to examine the folding pattern that 
is consistent with self-contact. Initially, by continuity, it 
will remain the ground state. However, with increased 
crowding one begins to force up the average curvature: 
above cpe = 8.27 the 3-fold possesses lower energy than 
the touching symmetrical 2-fold and the e-cone may flip 
from n = 2 to n = 3. Equivalent behavior is expected if 
(Pe is increased. To analyze the stability of our solutions, 
it is necessary to examine the second variation 5^H of the 
total energy functional. This is complicated by the fact 
that the local constraint of isometry has to be imposed on 
the deformations about the conical background. What is 
remarkable is that the calculation is tractable. One can 
show that 5^H is of the form f ds (j)C(j) where (j) = n- SX 
denotes the normal deformation of the surface. The relic 
of the isometry constraint is that f ds tvcj) = 0. The op- 
erator C = df^ ^dsVids + 1{V2 + ^F/0 is self-adjoint 
and of fourth order in ds. The potentials Vi and V2 are 
functions of k, and its derivatives. Using a decomposition 
of (j) into Fourier modes, one can determine the eigenval- 
ues of C for arbitrary cpe and n. They are all positive; 
e-cones free of self-contacts are stable. 
Conclusions. We have described the equilibrium states 
of a cone exhibiting a surplus angle. If we suppose that 
growth is slow compared to any viscoelastic timescale, 
the surface finds its equilibrium and the approach we 
have presented describes the evolution of the shape of a 
growing conical tissue. If the circumferential arc- length 
increases linearly with the geodesic radius, the surplus 



angle will remain constant; the cone will scale as it grows. 
Another mode of growth involves an increasing surplus 
angle. An initially flat disc will develop into a 2-fold, al- 
though a fluctuation may favor some higher energy state 
which we have seen is stable. However, if at some point, 
the surplus angle reaches (p^^^^{2) the surface will make 
contact which is costly energetically. At some higher 
value one would expect the surface to flip spontaneously 
into a 3-fold. This will continue through a 4-fold and 
so fifth with increasingly higher speed. Above V^e,max5 
however, self-contact becomes unavoidable. Internal lo- 
cal pressure will build up as the spherical volume occu- 
pied by the cone is packed more and more densely. 

What we have learned about the e-cone lays the founda- 
tion for understanding more general morphologies. If a 
disc surrounding the apex is removed, the cone can re- 
lax into some other flat geometry. Indeed one can easily 
verify with a paper model that the n = 2 ground state 
is unstable with respect to such deformations. These 
truncated cones can also be glued together to model sur- 
faces which are not flat: a surface of constant negative 
Gaussian curvature can be approximated by a telescope 
formed by such annuli. 
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